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A discrete quantum phase space formalism is used to discuss some basic aspects of the spin 
tunneling occurring in Fe8 magnetic clusters by means of Wigner functions as well as Husimi dis- 
tributions. Those functions were obtained for sharp angle states and symmetric combinations of 
QO . the lowest energy doublet - the particular tunneling energy doublet - with the application of an 

external magnetic field. The time evolution of those functions carried out numerically allows one to 
extract valuable information about the dynamics of the states under consideration, in particular, the 
CNj , coherent oscillations associated with the particular energy doublet. It also shown that an entropy 

functional can be constructed out of the discrete Husimi distribution and that the analysis of those 
coherent oscilations process allows one to clarify the meaning of the entropy in this case. 
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PACS numbers: 03.65.Ca, 03.65.Xp, 75.45.+j 



I. INTRODUCTION 



The basic ideas governing the discrete phase space picture of quantum mechanics of physical systems described 
by finite-dimensional state spaces have been long established, and the concept of discrete Wigner functions and 
<— j , the discrete analogue of the Weyl- Wigner transformations are now recognized as important tools in dealing with 
those systems [H, 0, Within this description, the quasiprobability distribution functions are writen in terms of 
a pair of discrete variables - for each degree of freedom - whose domains are finite sets of integer numbers which 
constitute a finite lattice that stands for a quantum discrete phase space. In this context, some attempts were also 
made in order to present a unified approach to the discrete quasiprobability distribution functions in the sense of 
CN) ■ obtaining the Glauber-Sudarshan, Wigner and Husimi functions as particular cases of s-parametrized discrete phase 
space functions [J, la, Q . In particular it was then shown that a hierarchical order among them through a smoothing 
process characterized by a discrete phase space function that closely resembles the role of the Gaussian function in 
the continuous phase space can be established Q . From an operational point of view, it is clear that if we have a 
discrete Wigner function, the corresponding discrete Husimi distribution can also be directly obtained. Furthermore, 
discrete quasiprobability distribution functions in finite phase spaces have potential applications in quantum state 



discrete quasiprobability distribution Functions m unite phase spaces nave potential applications m quantum state 
tomography [g, 0, quantum teleportation [|| [^, 

M M dLphase space representation of quantum computers 
[HI], open quantum systems [3], quantum information theory [laj. and quantum computation [l6[, where they are a 
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i natural tool to deal with the essential features of the inherent kinematics of the physical systems. 

In the present paper we intend to show that the discrete phase space approach presented before 0, 0, [(| [l7| can also 
be applied to the case of magnetic molecules, in particular to the Fe8 magnetic cluster [T^], where certain important 
quantum phenomena, such as spin tunneling, are currently considered because they are based on strong evidences 
coming from results of measurements with great resolution. That magnetic cluster has been deeply investigated in the 
last years [H, HO, Hl| and has been characterized as an useful tool in studying the importance of quantum effects in spin 
tunneling. Among the interesting properties it presents, it must be remarked that it has a J = 10 spin associated with 
its magnetic moment, and that the anisotropy constants characterizing its structure arc experimentally determined 
which allows one to construct a bona fide phenomenological Hamiltonian from which spectra and wave functions can 
be obtained. Besides, since it presents a well determined cross-over temperature below which the tunneling of the 
magnetic moment is dominated by the quantum effects [20j |. the experimental data obtained in that temperature 
domain can give inestimable information about that quantum phenomenon. The application of external magnetic 
fields on the Fe8 cluster also gives additional information about its quantum behavior and reveals interesting features 
of the energy spectra [22j ■ In this form, the quantum kinematical content as well as the proposed phenomenological 
Hamiltonian associated with this magnetic cluster provide the basic elements from which we can develop and apply 
the discrete phase space approach. As a consequence, we show that Wigner functions, as well as Husimi distributions, 
can be obtained that allow us to describe the behavior of the tunneling energy doublet and to study the effects of the 
application of external magnetic fields. Furthermore, after a Wehrl's type entropy functional has been constructed 
out of the Husimi distribution, a clear interpretation of that entropy is achieved through the study of the behavior of 
the tunneling energy doublet. It is shown that the discrete phase space approach results are in good agreement with 
the direct numerical calculations. 
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This paper is organized as follows. In Section II we present a brief review of the discrete phase space approach, 
including the procedures through which we obtain the discrete Wigner functions as well as the Husimi distributions. 
The procedure of getting the time evolution of the Wigner function is also presented. The application of the formalism 
to the Fe8 magnetic cluster is carried out in Section III, and Section IV is devoted to the conclusions. 



II. DISCRETE PHASE SPACE FORMALISM: A BRIEF REVIEW 



The possibility of describing nonclassical states and expectation values of quantum systems through a quantum 
phase space picture has been widely explored in the literature mainly in what concerns those physical systems described 
by continuous variables [23| . In this section we want to recall the main aspects of a possible quantum phase space 
picture espccifically constructed for treating physical systems described by finite-dimensional state space. 



A. Wigner functions and the Liouvillian 

In the past, we have proposed a mapping scheme that allows to represent any operator acting on a finite-dimensional 
state space as a function of integer arguments 0, [ijj on a discrete iV 2 -dimensional phase space; in other words, we 
have extended the well-known Weyl- Wigner procedure [23| of operators mapping for finite iV-dimensional Hilbert 
state space cases. In this form, we are then able to study the quantum properties of a given system of interest through 
the mapping of the relevant operators and studying the corresponding mapped functions. Of course we can also 
obtain discrete Wigner functions that represent state operators upon a discrete 7V 2 -dimensional phase space Q. In 
this connection, it is direct to see that we can also construct the Husimi distribution function, as well as the analogue 
of the Glauber-Sudarshan function if this is desired [1, Q . 

From a general point of view, let us start by considering that, for the systems under consideration, the discrete 
Weyl- Wigner transformation of a given operator O into its corresponding representative in the discrete quantum 
phase space is given by a trace operation 0, H3] 

O (m, n) = ~j^Tr [G^ (m, n) O] . 

In this expression G (rn, n) is an operator basis which guarantees the mod N invariance in the discrete phase space, 
and is given by 

G(m,n)= J2 ^^eM fc < w -^( mfe +"')], (1) 



where S (k, I) is the symmetrized, orthonormal and complete operator basis introduced by Schwinger [24[ and — (N — 
l)/2<£<(N-T)/2. 

It is also well known that the time evolution of a quantum system can be described by the von Neumann-Liouville 
equation for the density operator. Now, with the help of the expression for the mapped commutator of two operators, 
it is straightforward to obtain the mapping of that dynamical equation on the corresponding discrete quantum phase 
space equation. As already shown before ^he von Neumann-Liouville time evolution equation (we will consider 

hereafter h = 1) 
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is mapped onto 
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,—W(u,v;t)= £{u,v,r,s)W{r,s;t), (2) 



where 



dt 

r,s=- 



£(u,v,r,s) = 2i ^ 



h(m, n) 



N 4 



X sin 



e [in<f>(a,b,c,d;N)] 



[jf V>c -ad) 

x e { 2 N i [a(u-m)+b(v-n)+c(u-r)+d(v-s)]} ^ 
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is identified with the discrete mapped expression of the Liouvillian of the system, h (m, n) stands for the mapped 
expression of the here assumed time-independent Hamiltonian of interest, and the phase $(a, b, c, d; N) guarantees the 
mod N invariance. The function of integers W (r, s; t) is the discrete Wigner function which is obtained by mapping 
the density operator defined in the finite-dimensional state space; for instance, for a pure state 

p(t)=\i,{t)){i>{t)\ 

we have 

W (m, n; t) = ±Tr [G* (m, n) \ $ (t)) (V (t) |] . 

With these assumptions and considering the time evolved density operator to be written in the form of a series [26| 

p(t) = p(t )-i{t-t )[H,p{t )] 



2 (t-t ) 



2 



2! 

we see that its mapped expression is given by 



[H,[H,p(i )]]-..., (4) 



Z. 

W(u,v;t) = J2 W(^o){C ] C ] 



+ (-i)(t-t )£(u,v,r,s)+ { i)2( * *° )2 

x ^ a;, y) £ (x, y, r, s) + ...}, (5) 

where now W (u,v;t) is the solution at time t of the von Neumann-Liouville equation for the Wigner function. In 
expression (O those terms between curly brackets can be identified as the mapped propagator that associates the 
initial Wigner function at time to with a time evolved Wigner function at time t. We also see that the time evolution 
expression is constituted of sums of products of the array associated with the Liouvillian and the discrete Wigner 
functions taken over the discrete phase space. Then, from a numerical point of view, to calculate the time evolution 
of a particular density operator we only have to construct the arrays associated with the Liouvillian and the Wigner 
function of interest at time to- 

B. The Husimi distribution function 

From another perspective, in what concerns the mapped expression of the density operator, we can also have at our 
disposal a discrete Husimi function. The formal way of obtaining this function in the present context was discussed 
in previous papers in which the continuous Cahil-Glauber formalism (27j was extended for finite-dimensional state 
spaces 0, H| . Let us briefly review how to obtain the discrete Husimi distribution from a previously given discrete 
Wigner function. This can be accomplished by means of an s-parameterized mod(A r ) operator basis, which is given 
as 

T^) (m n) = V c [«^i^;'V)-¥(im+er.)] s(s) faj) 

Here the new mod(iV) invariant operator basis, (m,n), is defined through a discrete Fourier transform of the 
extended operator basis, namely 

s« fa, o = OP s fo, 0, (6) 

where K (rj, £) is a bell shaped function expressed as a product of Jacobi theta functions 

K(fi,Q = {2[0 8 (Q|ia)0 3 (Q|4ia) 

+ ^ 4 (0|m)??2(0|4m)]} _1 |i? 3 (7rar?|ia)i? 3 (7r<|m) 

+ tf 3 (7ra77|ia)??4(7ra£|ia)e ( "" )) 
+ i? 4 (7ra7 ? |m)?93(7ra^m)e ( ^ ) 
+ d i (7Ta V \ia)M^\ia-)e [l7 ' {T ' +i+N)] } 
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with a = (2Ny L and S (??,£) are the Schwinger basis elements already mentioned before. The complex parameter s 
that appears in the above expression is limited to the interval |s| < 1. 

In analogy with the decomposition expression of an operator into an operator basis, what involves its representative 
defined on the discrete quantum phase space, we have now 

1 1 

= — N (m,fi)TW(m,n). 

m.n——£ 

In particular, if we choose O = p we see that for s = this expression gives the mapping leading to the Wigner 
function, while for s = — 1 it leads to the Husimi distribution. 

Within this formalism we also have the expression that allows us to obtain the Husimi distribution directly from 
the Wigner function, namely 

t 

H{m,ri)= ^2 E ( m ' ~ m > n ' ~ n ) w ( TO '' n ') ' 

m' ,n' ——£ 

where E (m! — m,n' — n) characterizes the mapping kernel of the transformation, and it is given by 

E(m! -m,ri -n) =Tr T (0) (m, n) T^ 1 ' (to', n') . 

In what concerns the time evolution, it is then clear that, once the Wigner function is given at any instant of 
time, the Husimi distribution can be obtained straightforwardly in the same way. Furthermore, in what relates to 
the system state, since the Husimi distribution has a well-behaved shape, we expect to obtain a clearer picture of the 
physical process of interest through its use. 



III. APLICATION: FE8 MAGNETIC CLUSTER 



The discovery of magnetic molecules, in particular the Fe8 magnetic cluster, which was first synthesized by 
Wieghardt et al fl8| . was a great contribution for the study of quantum effects in the dynamics of spin tunneling 
[28| . Under proper temperature conditions, the pure quantum effects in the dynamics of the magnetization become 
manifest in some measurements performed with this molecule, specially due to the very well measured constants 
characterizing its structural properties. The study of the orientation of the magnetic moment with these clusters also 
benefits from the fact that the presence of an external magnetic field may enhance some quantum transitions that 
allow one to test basic hypothesis about the quantum spin tunneling. 

The main experimentally determined results related to the Fe8 magnetic cluster indicate that it is a system with 
a J = 10 spin - that implies in a 21-dimensional state space for this degree of freedom- with an observed barrier of 
about 24 K . Furthermore, it was verified that below 0.35 K the relaxation of the magnetization becomes temperature 
independent [20l |. and that when an external magnetic field with intensity that is given by integer multiples of 
Aff|| ~ 0.22 T is applied along the easy-axis of the cluster (z-axis) the hysteresis curve has a behavior ^lj which can 
be interpreted as if there is an energy matching of states on both sides of a double-minima barrier, similar to what 
had already been proposed in another magnetic molecule, namely the Mnl2ac cluster [29| . From the experimental 
data the proposed Hamiltonian for the Fe8 cluster (sometimes called giant spin model) in the presence of external 
magnetic fields has the form 

H = Ho + H!, (7) 

where the term 

H =£>J* + | (J* +J 2 _) 

constitutes the part that is established by the measured properties of the molecular structure, and the remaining ones 

H 1= AJ Z + B (J + + J_) + C (J+ - J_) 

give the contribution that is associated with the external magnetic field. The anisotropy constants are plj D/kB = 
—0.275 K and E/ks — 0.046 K, where ks is the Boltzmann constant, while A = g^igH^, B — gyuB cos aH± , and 
C = g[iB sinaiJj_ are the parameters associated with the magnetic field intensity along the z-axis (easy-axis), y-axis 
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(medium axis) and x-axis (hard-axis) respectively. Furthermore, this model Hamiltonian is expressed in terms of the 
angular momentum operators that obey the standard commutation relations 

[J±,J Z ]=TJ±, [J+,J-] = 2J z , 

and its numerical diagonalization can be directly carried out within the 21-dimensional set of \j, to) states, so that 
those energy eigenvalues thus obtained can be useful as reference values. 

Now, the main requirement for the use of the proposed description of physical systems with finite-dimensional 
state space consists in characterizing the starting space state such that the quantum phase space mapping scheme 
is completely established. In the present case we have to establish a connection between our approach and the 
mathematical description of the proposed model for the Fe8 magnetic cluster. This connection, which is then the 
key element in order to allow the use of the quantum phase space description for the Fe8 cluster, is obtained just by 
choosing 

\u k ) = \j, k) 

with the associated eigenvalue equation following the prescription proposed by Schwinger [2~ij 

U K) = exp ^TriAj \j,k), (8) 

where U is an unitary operator and \uk) its eigenstates. Furthermore, 

3z\j,k) = k\j,k), 
J 2 |j,fc) = j(j + l)\j,k), 



J±|j,fc) = VtiTk)(j±k + l)\j,k±l), 

as usual, where — 10 < k < 10. 

Therefore, we have a 21-dimensional space of the U eigenvectors as given above, and consequently we also have 
another 21-dimensional space, of eigenvectors of V, the complementary unitary operator also proposed by 

Schwinger. It is worth noting that the set of states {|t>z)}, where — 10 < I < 10, is in fact a basis in this finite space 
since it is orthonormal and complete, thus it is also a bona fide basis. Since it was also shown [2~ij that the set of 
states {|fz)} is obtained through a discrete Fourier transform of the basis states {|itfc)}, the physical interpretation of 
the states is then direct: they are associated with the polar angle orientation of Fe8 cluster spin. 

Based on this association, we can now use the discrete Weyl-Wigner formalism discussed before to describe the 
quantum behavior of the Fe8 cluster in a discrete 21 2 -dimensional phase space. Thus far the kinematical content of the 
description has been already set up, then we can go one step further and obtain the mapped expression associated with 
the Fe8 cluster Hamiltonian ([7]). Following the mapping scheme presented before the discrete phase space expression 
for that Hamiltonian is 

h (to, n) = Am + Dm + EQ (to, n) cos 

10 e ^i[ r ( m _ fe _l)_ n ] 



E 



2N 

k,r=-10 



x V U — k)(j + k + 1) (B cosa + C sin a) 

+ E - 



fc,r=-10 



2N 



x y/(J + k) (j : - k + 1) (B cos a - C sin a) , 

where 

Q (to, n) = \J (j + m) (j — to + 1) (j — to) (j + m + 1). 



With this result, it is straightforward now to find the expression for the Liouvillian. After a direct but tedious 
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calculation, the discrete phase space expression for the Liouvillian is obtained, and it is written as 

2 1 



C(u,v,r,s) = ^2 (J-(Am+ Bm 2 ) sin (jjad) 

m,a,c,d 
iTr&(a,0,c,d:N) 

e 

Cy/(j + m)(j -m + 1) sin \j-(c - ad)] 

X eW^+^g 2 ^- c(r-a) e iwS>(a,l,c,d;N) 

- IV C? - m ) U + m + l) sin [iL(c + ad)] 

x e N 1 'e N e v ' ' ' ' 



x e 

+ 



+ 



Q(m, n) 



:,d;N) 



x |sm (2c — ad)Je N e ^ ' ' ' 
— sin ^— (2c + aa)Je e v ' jj 



Jp- [a(« — m) + c(u — r) + d( 



(9) 



Having obtained the discrete Wigner function and using the mapped series, Eq. ((5|), we are then able to compute 
the temporal evolution of a state given at r = with the expression above, Eq. ©. 

In what concerns the Wigner function at r = 0, we may choose those that may be relevant for our purposes. At 
first, we start, for simplicity, from the density operator p for a pure state 

where can be expressed in terms of the eigenstates of the operator U, Eq. ([5]), which are also, due to the adopted 
formalism, eigenstates of the angular momentum operators J 2 and 3 Zl namely 

10 

W)= E <w>- 

l=-10 

As has already been shown elsewhere [30| . the Wigner function in this case is written as 

10 



W(m,n) 



1 

N 



E exp 



r,s=-10 



2iri 



(mr + ns) 



9 (r, s) , 



(10) 



where 



10 



{k+s} exp 



fc=-10 



2ni /, s\ 



and the symbol {A; + s} stands for the sum cyclically restricted to the interval of labels, e.g., /c max + 1 = fcmin- As 
can be directly recognized, a great advantage of using the Wigner functions in this form is that, by diagonalizing the 
Hamiltonian in the angular momentum state basis {\j, k)}, we can obtain the coefficients Ck to construct the Wigner 
function through Eq. (jlOjl . 



A. Sharp angle state 

A first illustrative case that can be directly treated is a sharp angle state associated with the physical configuration 
where, at r = 0, the magnetic moment of the Fe8 cluster points to a particular angle {8,(p}. The parallel and 
transverse magnetic fields are, in particular, Hu = 0.01 T, and H±_ = 2.0 T respectively, being H j_ oriented at an 
angle of tp — 7r/4 relative to the x-axis, and the initial state, a sharp angle, is peaked at n — 0, where 6 n = j^n with 
— 7r < 6 n < 7r, and — 10 < n,m < 10. The correspondent sequence of the time evolved Wigner functions is depicted in 
Fig. 1. 
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(a) 0>) (c) 




FIG. 1: Time evolved Wigner function calculated with time steps of At = 0.05 x 2ir K" 1 and external magnetic fields H\\ = 0.01 
T and H± — 2.0 T oriented at an angle of <j> = tt/4 with the a>axis. The labels m and n correspond to the finite-discrete angular 
momentum and angle variables respectively displaced 11 units of their original intervals. Figure (a) corresponds to the Wigner 
function at r = 0. Figures (b), (c), and (d) depict the Wigner functions for r = 0.4 x 2n A' -1 , 0.8 x 2n K" 1 and 1.2 x 27r K' 1 
respectively. 

Now, since we have the time evolved Wigner functions for the peaked angle state, the associated Husimi distributions 
can be directly obtained, and are shown in Fig. 2 

In general, if the Wigner function evolves in time not clearly showing the behavior of the physical system, the 
Husimi distribution, on the other hand, can afford for a better visualization of how the state evolves over the discrete 
phase space. In particular, it is to be noted the smeared behavior that the Husimi distribution presents. 

Also we see that, after a transient time - during which the initial distribution loses its identity - the peak of the 
distribution passes recurrently on the same regions of the 9 = [0, ir] sector of the discrete phase space because of the 
quantum correlations and the presence of the transverse magnetic field in the xy-plane. 

B. Symmetric combination 

As a second case we consider the symmetric combination, namely, 




where the labels i an j characterize the -Fe8 cluster energy eigenstates of which the combination is made of. Here we 
will be interested in the energy doublet constituted of the ground state and its neighbor partner. This choice is based 
on the fact that this doublet is the proposed candidate for showing the spin tunneling properties when the Fe8 cluster 
is submitted to temperatures below the well established crossover temperature. Furthermore, it is also known that 
the energy gap between them is widened when an external magnetic field - parallel to the easy-axis - is applied on 
the cluster. In particular, this study is of great interest when this external magnetic field, H\\, is an integer multiple 
of 0.22 T because of the matching of energy levels as pointed out before. 

For those states it is also simple to obtain the Wigner function by means of Eq. (fTTj)) . namely W s (m,n), and 
propagate it in time. In Fig. 3 we depict a sequence of time evolved Wigner functions for H\\ = 0.11 T introduced in 
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(a) (b) (c) 




FIG. 2: Time evolved Husimi function under the same conditions indicated in Fig. 1. 



order to induce a convenient widening of the energy gap, and H±_ = 0.0 T. For this situation we have, through the 
diagonalization of the Hamiltonian, the reference values E g = —29.01745 K and E\ — —26.06441 K which are the 
ground and the first excited states respectively. The energy gap beetwen those states is AE re f = 2.95304 K. 

First of all we observe that the staggered behavior of the Wigner function, mainly visible in the central region of the 
discrete phase space, is due to parity effects related to the angular momentum variable. Moreover, what is important 
here is that the initial configuration is totally recovered after r « 2.15 x 2ir A"" 1 , as can be promptly recognized; this 
periodicity with time of the Wigner function indicates that we have a coherent oscillation that is associated with the 
magnetic moment quantum tunneling. 

If, instead of the Wigner function, we have the family of the corresponding Husimi distributions, as shown in Fig. 
4, we immediately see that the staggering was smeared. Also, the oscillation can be clearly seen to occur along the 
angle sector of the phase space such that the periodicity is obtained and agrees with that observed with the Wigner 
function, in other words, the smearing process does not alter the oscillatory character of the motion, so that the 
periodicity is the same. 

Note that, the periodicity of the Wigner function associated with that combination of states could also be obtained 
from the time correlation function defined as 

P if ( r ) = (™> n ; °) W f {m, n; r) . 

m,n 

A calculation performed with the same external magnetic field as before (H\\ = 0.11 T and H± = 0.0 T) gives Pif (r) 
as shown in Fig. 5. Once again we obtain for the symmetric state a periodicity in time of r w 2.15 x 2ir K . Of 
course the oscillation period can be obtained from Fig. 4 as well from Fig. 5; both coincide. 

Now, as an interesting result we see that we can associate this oscillation period, t (H\\), with the energy gap 
between the ground state and its neighbor partner - which give the first energy doublet - by just using 

where we can associate a frequency ui (#||) with the oscillation period, and AE is the energy gap given as a 
function of the magnetic field intensity. It is important to emphasize that this relation is dependent on the magnetic 
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FIG. 3: Time evolution of the Wigner function generated by a symmetric combination of the lowest energy levels. The time 
step was taken as At = 0.05 x 2n K~ , and the external magnetic fields intensities were Hn = O.lf T and H± = 0.0 T. Here 
the labels m and n also are the discrete angular momenta and angle variables. Figure (a) corresponds to the Wigner function 
at r = 0. Figures (b), (c), and (d) depict the Wigner functions for r = 0.75 x 2tt K~ x , 1.55 x 2ty K" 1 and 2.1 x 2n K~ 1 
respectively. 




FIG. 4: Time evolved Husimi distribution under the same conditions indicated in Fig. 3. 
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T [27TK- 1 ] 

FIG. 5: Time correlation function calculated with time step of At = 0.05 x 2n K . The external magnetic field intensities 
are H\\ = 0.11 T and H± = 0.0 T. 

field intensity, and that here we want to consider only week magnetic field intensities effects on the tunneling spin 
energy doublet. For the values of the above considered external magnetic field and using r = 2.15 x 2n K -1 , we 
obtain from Eq. (fTTj) AE = 2.92241 K, while the reference value is AE re f — 2.95304 K, thus showing a deviation of 
1.04%, which indicates a good agreement. 

Still along the same line we have just done, and in order to further clarify the tunneling process ocurring in this case, 
we can also use an angle-based approach developed in recent years (3lL [32j to describe this same physical situation. 
In that approach we obtain an effective Hamiltonian, written in terms of a potential and an effective mass function, 
that can describe the present situation and that can also complement the present discussion. To highlight the main 
spin tunneling characteristics visualized in Figs. 4 and 5 we plot the associated potential function, Fig. 6, 

V{9) = (D + E)J(J + l)cos 2 (<9) 

- g( i B H\\y/J(J + 1) cos(fl) - EJ{J + 1), 




-10 1 2 3 4 
6 [rad] 



FIG. 6: Potential function of the Fe8 cluster obtained through an angle-based approach. The magnetic field intensities are 
ffy = 0.11 T and H ± = 0.0 T. 

where the Fe8 parameters are those already presented before. We then see that the peaks of the wave functions which 
generate the Husimi distribution - associated with the symmetric combination at r = - occur at both deep minima 
at 9 = and 6 = tt of the potential. It is then immediate to verify that the leakage of the Husimi distribution into 
angular regions that would not be classically accessible indicates that tunneling is ocurring. 

C. Entropy Functionals 

As is well known, the Husimi distribution is an important tool to infer the quantum information associated with 
the particular quantum state of the physical system it represents. In this connection we can in particular write a 
Wehrl's type entropy expression 33] in terms of the Husimi distribution 



S(H) = --Tr[H log [H]] 
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Now, since we can consider the Husimi distribution as a matrix over the discrete phase space - it assumes a positive 
value on every site of the phase space -, we can write the entropy as [34j 

1 1 

S(H;t) = -- Y, H Kn;r) log (W(m,n;r)), (12) 

m,n— — £ 

where the sums are taken over the entire phase space, and it is immediate to see that it is non-negative. 

Taking that expression as our starting point, we can extract its inherent physical interpretation by just choosing a 
particular quantum state of which the essential features of the corresponding Husimi distribution have been previously 
analysed. Moreover, since we are interested here in Husimi distributions associated with simple pure states, we will 
focus our attention on the time evolution of the entropy functional since it can give us essential information about 
the time evolution of the particular state all over the discrete phase space, in particular the indications of coherent 
oscillations associated with spin tunneling. 

As an emblematic case to be treated here, we consider the symmetric combination of the lowest energy doublet of 
the Fe8 cluster submitted to an external magnetic field parallel to the easy-axis. As has been seen before, in the time 
evolution of this state a period of oscillation is recognized that characterizes the energy gap in the spectrum. On the 
other hand, this same fact can also be studied from the entropy point of view, what will give us some indications of its 
interpretation [36]. Performing the numerical calculations for ify — 0.11 T and H± = 0.0 T, we get the corresponding 
value of the entropy at each step of the time propagation. The oscillatory character of the entropy thus obtained is 
depicted in Fig. 7. 

0.217 
0.2165 

"p 

£ 0.216 
0.2155 



0.5 1 1.5 2 

t[27tK- 1 ] 

FIG. 7: Time dependence of the entropy function of the symmetric combination of the two lowest states using the same 
parameters as in Figs. 4 and 5. 

Based on these results we verify that the Wehrl's type entropy in this case measures how the Husimi distribution 
spreads over the angle and angular momentum sectors of the discrete phase space. But we have to consider that the 
entropy is a minimum whenever the Husimi distribution is maximally concentrated in a particular region of the phase 
space, and is a maximum when the distribution spreads over the phase space. Since the entropy does not distinguish 
negative or positive sectors of the phase space, we have to consider for the period of oscillation the time it takes for the 
system to start from a minimum in a particular region of the phase space and to return to that same region. Taking 
this into account, the energy gap can be calculated from the period thus obtained and agrees again with the reference 
value within an error of 1.95%. Thus, by observing the oscillation of the entropy associated with the symmetric - 
or antisymmetric - combination of an energy doublet we can infer the gap energy information associated with the 
magnetic moment tunneling. 

Besides, we can also use the mutual correlation functional [H, [3f| as well 

1(H; t) := S(J; t) + 5(6; t) - S(H; t) > 0, (13) 

where S(TC;t) corresponds to time-dependent entropy defined in Eq. (fT2"|) . S{J]t) and S(Q\t) represent the marginal 
entropies which are related to the respective marginal distributions J{m\i) and <d(n;t) [q. |36|] . obtained by taking 
the trace of the corresponding Husimi distribution. With this functional we recover exactly the same behavior for the 
symmetric combination as that presented in Fig. 7, which allows us to infer possible correlations between the states 
and also calculate the energy gap as it was done through the Wehrl's type entropy. 

IV. CONCLUSIONS 

Based on a formal approach aiming at the construction of a discrete phase space picture of quantum mechanics 
related to finite-dimensional state spaces developed before, and briefly reviewed here, we addressed the problem of 
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spin tunneling in the Fe8 magnetic cluster. Our formal approach was shown to hold also in this particular case of 
finite-dimensional quantum system, and starting from the phenomenological Hamiltonian describing the magnetic 
cluster with or without an external magnetic field we obtained the discrete phase space representatives of the relevant 
elements which govern the complete description of the system and its time evolution. The Liouvillian associated with 
the Hamiltonian - which propagates the density operator in time in the discrete phase space - was explicitly written. 
In what concerns the representation of the state of the Fe8 magnetic cluster we have shown how a discrete Wigner 
function, as well as a discrete Husimi distribution, can be directly obtained. 

The time evolution of those functions for the Fe8 magnetic cluster was studied in two different cases: a) sharp angle 
state with an external magnetic field, and b) symmetric combination of two energy eigenstates of the phenomenological 
Hamiltonian; in particular, we chose the lowest energy doublet. In both cases the Wigner function, as well as the 
Husimi distribution, were calculated and relevant information concerning the dynamics of the spin was extracted 
from them. With respect to the time evolution of the symmetric combination of states, it was also shown that, by 
analysing the periodicity of its oscillatory behavior, a period could be found that allowed us to calculate the energy 
gap associated with the tunneling doublet of states from which we started. 

Finally, we introduced a Wchrl's type entropy functional based on the discrete Husimi distribution that allowed us 
to obtain a measure of the information related to how the Husimi distribution is distributed over the discrete phase 
space. Well concentrated distributions in phase space are associated with small entropies, while the opposite produces 
greater entropies. The coherent oscillations verified in the symmetric combination of the lowest energy eigenstates 
thus reflect an oscillation in the associated time dependent entropy. Therefore, it was also possible to extract the 
value of the energy gap of the energy doublet through an analysis of the entropy. Furthermore, the values obtained 
via the mutual correlation functional lead us to the same physical interpretation as that obtained with the Wehrl's 
type entropy functional. 

The results presented here corroborate all the basic formal scheme proposed as well as the procedures adopted here 
to describe spin systems, in particular the Fe8 magnetic cluster, in the corresponding discrete phase space. 
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